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Abstract 

Supersymmetric SO(IO) grand unified models with renormalizable Yukawa couplings in- 
'sj" ' volving 10, 126 and 120 Higgs fields have been shown to give a very economical theory for 
[P^^ ■ understanding quark-lepton flavor in a unified framework. In previous papers, we showed how 
CN ! nucleon decay can be suppressed in these models without invoking cancellation, by choice of 
^ I Yukawa flavor texture within a type II seesaw framework for neutrinos that explains all mixings 

^ I and masses including the recently observed "large" ^13. In this follow-up paper, we extend our 
earlier work to the case of type I seesaw and show that the recently measured "large" ^13 can 
be accommodated in this case while suppressing proton decay. We then point out that the two 
cases (type I and II) lead to different testable predictions for B{fi — e + 7) and B{t /i(e) +7) 
as well as different flavor final states in nucleon decay. In particular, we find that for the type I 
seesaw case, B{t — )■ yU + 7) can be observable while at the same time suppressing — )■ e-l-7), 
whereas in the type II seesaw case, i?(r —)■ /i + 7) is always suppressed whereas — )■ e -|- 7) 
is observable. 
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1 Introduction 



The last decade has seen incredible experimental progress in the field of neutrino physics. Since 
the discovery of the oscillations of atmospheric neutrinos in 1998 followed by the confirmation 
of the solar neutrino oscillations discovered by Ray Davis in mid-sixties, many parameters 
describing the neutrino masses and mixings have been measured in various solar, atmospheric 
and accelerator and reactor experiments. The most recent such measurement is that of the 
remaining mixing parameter 613 by reactor as well as accelerator experiments [H [21 [3l HI |5] 
which has created a lot of excitement in the theory community due to its implications for 
physics behind the neutrino masses. Attention has now turned to measuring the other missing 
pieces of informations: Dirac versus Majorana nature of neutrinos, the CP violating phase S 
in the lepton sector as well as the mass hierarchy of neutrinos. The global analysis of the 
parameter fit [0] seem to have some preference for CP phase ~ vr although it is far too early to 
take this seriously. The "large" value of the 13 mixing angle (around sin^ia ~ 0.15) however, 
has raised the hope that the CP phase can be measured in near future. 

What have we learnt about the physics behind the neutrino masses from these observations? 
The first challenge is to understand why the neutrino masses are so small and the second is to 
see if the same framework that explains the small neutrino mass also simultaneously explains 
the observed mixing pattern. A popular paradigm for this seems to be the so-called seesaw 
mechanism [7j where one extends the standard model by adding heavy Majorana right-handed 
neutrinos (known as the type I seesaw) or heavy SM Higgs triplets (known as the type II 
seesaw) Generic models based on seesaw are quite successful in achieving the first goal 
i.e. small neutrino masses. It is however more challenging to get an understanding of the 
mixings without further assumptions. The two approaches which have been used extensively 
are: (i) leptonic symmetries and (ii) grand unified theories. The closeness of the leptonic 
mixing parameters to group theory-like numbers (e.g. ^) have been the major driving 
force behind the symmetry approach. Most models with symmetries however predicted the 
third mixing angle 613 = 0, which have now been disproved by experiments as noted. While 
there are still a number of leptonic symmetry scenarios which can lead to large value of 613 [9] , 
the fact that 



sin 613 ~ A — — 1^ ~ Ocabibbo (1) 
\ Am23 

may be suggestive of a "big" picture that unifies quarks and leptons such as those based on grand 
unified theories. In particular, if normal neutrino hierarchy is established, that would imply a 
similarity between the quark and the lepton sector, that would be consistent with expectations 
from grand unified theories (GUTs) and would further motivate the GUT approach to all 
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flavors. 

One class of grand unified theories that have been extensively explored and appear quite 
successful in providing a unified description of quark-lepton flavor i.e. giving observed values of 
the neutrino mixings and masses, while being consistent with observed quark masses and CKM 
mixings is the one based on the group SO(IO) with or without supersymmetry and using 126 
fields to break the B — L subgroup of SO(IO). This class of models uses only renormalizable 
Yukawa interactions involving a 10 and 126 [TUl E] with/without an additional 120 [121 [13] 
to explain both the quark and lepton flavor puzzle. The renomalizability property restricts the 
fermion Yukawa couplings to have only a small number of parameters (without the need for 
adhoc symmetries) describing both the quark and the lepton sectors simultaneously so that 
the models become quite predictive. Indeed these models led to the prediction not only of the 
relation between neutrino mass ratios and Cabibbo angle in Eq.([T]) but also the value of 6'i3 
several years prior to its measurement |14j without the need for 120 fleld. 

One of the key predictions of the supersymmetric (SUSY) GUT models is the enhancement 
of proton decay rate due to the presence of new supersymmetric contributions from color 
triplet Higgsino exchange. These contributions are known to severely constrain the nature of 
the these theories and clearly if the above SO(IO) neutrino models are to be taken seriously, 
their consistency with current lower limits on proton lifetimes must be examined for both 
the decay modes p — )■ e^vr^ where the key input is the value of the uniflcation scale and the 
characteristic SUSY mode p — )■ K~^u, where the key input is the nature Higgs fields that give 
mass to fermions. For the case of minimal SO(IO) models with 10, 126, with type II seesaw 
for neutrino masses, this consistency requires cancellation between parameters describing the 
Higgsino mixings and masses [15j. Since there are many proton decay modes that have been 
constrained by experiments [16] and all parameters of the model except for three Higgsino 
mixngs are already determined, the fact that there is a consistent picture, is highly nontrivial. 
However this works only for low tan/3 regime of the MSSM parameter space. It was subsequently 
pointed out [171 HH] that once the 120 is included, the situation improved quite a bit within 
the type II seesaw framework i.e. by choice of an appropriate fiavor structure for the different 
Yukawa coupling matrices, one can not only give a simultaneous fit to fermion masses but also 
suppress proton decay to the desired level, for both small and large tan (3, without need to 
invoke cancellation. It was shown that in the model with 120, one can obtain ^13 as large as 
0.15, which is in accord with the recent measurements. 

An important question remained as to whether the idea of suppressing proton decay by 
choice of fiavor structure works in models with 120 when type I seesaw is picked for under- 
standing neutrino masse. The second question is: is their a way to distinguish between the 
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type I and type II seesaw models of this kind. Clearly, there are two obvious avenues to ex- 
plore in this connection. One can look at the flavor structure of the proton decay flnal states 
and secondly the predictions for rare lepton flavor violsting decays of type — > e + 7) and 
B{t —7- yu(e) + 7). We investigate both these questions in this paper. 

The reason to suspect that lepton flavor violation can probe different models is that in SUSY 
seesaw models all leptons are accompanied by their bosonic partners, the sleptons and their 
mixings can lead to large flavor violating effects at low energies. To prevent excessive flavor 
violating effects, one generally assumes that at some high scale, all slepton masses are equal. 
However as we extrapolate the theory down to the weak scale, slepton mixing are generated by 
the large neutrino mixings hidden in the Dirac Yukawa couplings of right-handed neutrinos |19j . 
In generic type I seesaw models without additional inputs, these rates can be suppressed by 
simply changing the seesaw scale since these mixings are proportional to products of Y^, matrix 
elements and if the seesaw scale is lowered, seesaw formulae for neutrino masses demand that 
Yi, become smaller in absolute magnitude thereby reducing the — )■ e + 7 etc. rates. The same 
thing also happens in SU(5) GUTs, where the size of the Dirac neutrino Yukawa coupling is a 
free parameter because the right-handed neutrino is a gauge singlet and no deflnite prediction 
is possible. 

However, in predictive GUT models such as the minimal renormalizable SO(IO), coupling 
uniflcation and fermion mass flts determine all Yukawa couplings of the theory and in particular, 
the Yukawa couplings of the RH neutrinos. Furthermore the values of the RH neutrino masses 
are also predicted making the lepton flavor violation predictions more deflnite. There are also 
separate contributions to LFV amplitudes coming from type II seesaw Yukawa couplings [20] , 
which are not generally considered. Another contribution to LFV in all GUT models comes from 
charged fermion-quark-color triplet Higgs {c'^u'^Ht) coupling |2T] which is independent of the 
seesaw mechanism; it is also predicted in these models since fermion mass flts also determine 
those couplings. Moreover, since the flavor structure of the Yukawa matrices is related to 
the lepton flavor violations (LFV), and the proton decay amplitudes within a SUSY SO(IO) 
framework, these two predictions are in principle connected. We investigate this question in this 
paper for a speciflc realistic SO(IO) model. We note that although some early studies of LFV in 
SO(IO) models have been carried out in Ref. [22], where assumptions are made about expected 
relations between neutrino Yukawa matrix Y^ and in realistic models there may be signiflcant 
deviations from these relations. We flnd this to be true in models we are considering where in 
addition to new forms for the 10 couplings, there are important 126 Yukawa contributions to 
flavor mixing. 

The new results of this paper are: (i) we present a detailed analysis of the necessary texture 
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to obtain "large" ^13 and suppressed proton decay for type I seesaw case and in particular 
discuss how the relation of 13 mixing and the ratio of the squared mass differences, Eq.(IT]) 
emerges: (ii) we show that in the case of type II seesaw mechanism, zero value of 13 neutrino 
mixing is disfavored if we want to suppress proton decay via flavor structure; (iii) we also show 
that the measured value of 13 neutrino mixing prefers 5pmns ~ ti" for the CP phase in the 
neutrino oscillation paramaters and (iv) we compare the predictions of LFV for both type I 
and II seesaw and show how the required flavor structures in type I and type II seesaw cases can 
provide different predictions for LFV decays: in particular, we find that for the type I seesaw 
case, i?(r — /i + 7) can be observable while at the same time time suppressing — )■ e + 7), 
whereas in the type II seesaw case, B{t — )■ /i + 7) is always suppressed. 

This paper is organized as follows: in section II, we discuss the Yukawa matrices in SO(IO) 
model; in section III, we discuss the suppression of nucleon decay amplitudes; in section IV, we 
discuss the LFV constraints; in section V, we discuss the flavor structures and predictions for 
type II seesaw scenario; in section VI we discuss the flavor structures and predictions for type 
I seesaw scenario; in section VII we discuss the predictions for lepton flavor violation in type I 
and II seesaw scenarios; in section VIII we discuss the predictions for nucleon decay in type I 
and II seesaw scenarios and we conclude in section IX. 



2 Yukawa matrices for fermions in renormalizable SO (10) 

The Yukawa terms in the superpotential of the renormalizable SO(IO) model involve the cou- 
plings of 16- dimensional matter spinors ipi with Higgs fields belonging to 10 (denoted by H) 
and A, and D, representing the 126, and 120 dimensional representations, respectively and is 
given by: 

Wy = ^hijAi'jH + A + h[.ij,ij,D. (2) 

This equation holds at the GUT scale. In order to write the effective Yukawa couplings below 
the GUT scale, we extract the effective Hu^ fields of MSSM which are linear combinations of 
the Higgs doublet fields not only in H, A, D fields but also in other fields e.g. A and 210 fields 
used to break the GUT symmetry while maintaining supersymmetry down to the weak scale. 
We assume that the SO(IO) symmetry is broken down to the standard model gauge symmetry 
by the vacuum expectation values of A(126) + A(126) and the 210 Higgs field. Our conclusions 
below are independent of this i.e. the GUT symmetry could have been broken down by other 
fields such 54, 45 etc as long as they do not contribute to fermion masses. We also assume 
that only one pair of the linear combinations of the Higgs doublets remains massless (or more 
precisely, weak scale Higgsino mass, which is much smaller than the GUT scale) to break the 
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electroweak symmetry. Using this, we can write the Yukawa matrices that give rise to the 
fermion masses as those given by the hnear combination of the original h, f and h' couphngs: 



Yu = h + r2f + rsh', (3) 

Yd = n{h + f + h'), (4) 

Fe = n{h-3f + Ceh'), (5) 

Yu = h-3r2f + c,h', (6) 



where and Ce, are the functions of Higgs mixings, and h, f and h' matrices are the original 
couplings multiplied by Higgs mixings. More details about these equations can be found in 
Ref.[n]. 

The 126 Higgs Yukawa coupling includes both left- and right-handed Majorana neutrino 
couplings: 

ipipA D UAl + z/z/Ar, (7) 

where denotes a S'U{2)l triplet in the 126 representation. As a result, in the case of 
triplet-part dominant type II seesaw neutrino mass [8], the neutrino mass is (approximately) 
proportional to the coupling matrix /. It is clear from the above equation that there is an 
intimate connection between the lepton and the quark sector since the same Yukawa coupling 
matrix / appears both in the charged fermion sector as well as the neutrino sector due to the 
seesaw mechanism. We assume CP conservation prior to symmetry breaking so that h and 
/ are real symmetric matrices and we can choose a basis so that h is diagonal. The matrix 
h' is imaginary and anti-symmetric so that the total number of coupling parameters in the 
theory (prior to any assumption about nucleon decay) is twelve. Combined with the other 
parameters for the model i.e. (adding in ri,Ce,u) and the vbl which gives the overall scale of 
neutrino masses) the total number of parameters is eighteen. This is smaller than the number 
of observables i.e. 13 in the charged fermion sector and six in the neutrino sector excluding 
the Majorana phases. So the model is predictive. When combined with the flavor ansatz for 
suppressing nucleon decay mentioned above, the number becomes fifteen and the predictive 
power increases as we saw in Ref . [17] . 

3 Suppression of nucleon decay amplitude 

In this section, we review the flavor suppression mechanism for nucleon decay in this model 
proposed in [17J. The dominant proton decay in SUSY GUTs arise from the exchange of color 
triplet higgsinos which are part of the 10, 126 and 120 fields. The effective dimension five 
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operator induced by them can be written as 



(8) 



where the coefficients Cl,r are functions of /i, / and h' couphngs: 



^ijkl 




lijkl 
R 



E -T^iXalh - XaJ + V2X,ih'),,{Y,,h - (F,5 - V2Yae)f + V2{Ya3 - Ya2)h')kl, (10) 



where Mt^ are masses of colored Higgs triplets, and X and Y are the colored Higgs mixings 
(Details are given in Ref.|17j). 

As is noted in the introduction in order to suppress the nucleon decay amplitudes, one has 
to consider 

1. Choice of SUSY breaking parameter (heavy squarks being preferable). 

2. The (lightest) colored Higgs triplet is heavy and/or the colored Higgs mixings are small|l| 

3. Special structure of Yukawa matrices. 

All the above requirements may work in tandem to satisfy the current experimental constraints 
on proton lifetime without invoking cancellation. We will concentrate on the item 3 in this 
paper, and consider what kind of flavor structure is favorable to suppress the nucleon decay 
operators, while at the same time giving correct prediction for the fermion masses and mixings. 

The important features that arise in the discussion of suppressing nucleon decay operators 
by a flavor structure are as follows [171 123] : 

1. In the current setup, there are multiple pairs of Higgs fields, and thus, there is freedom 
to cancel the decay amplitude. However, the required cancellations in both Cl and Cr 
operators are large, especially for large tan/3. Sometimes, Cr operators are ignored since 
they are suppressed by Higgsino dressing rather than gaugino dressing. It is true that the 
contribution from Cr to the nucleon decay amplitude is rather suppressed compared to 
Cl- However, in the case that the 1st generation masses are obtain by a choice of h and / 
coupling and hn, fu ~ ya (down quark Yukawa coupling), the Cr operators are far from 
being small. Indeed, even in the minimal-type of SU(5) model, when only one pair of 
Higgs fields couple to fermions and the dimension five operators are roughly in the form 

^We investigate what kind of SO(IO) breaking vacua is preferable to obtain heavy colored Higgs in the current 
scheme in Ref.|l8). 
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of Irfl^n, and leads to trouble due to the fact that both Cl and Cr cannot be suppressed 
simultaneously from the flavor structure [17] unless the colored Higgs spectrum is extended 
|18] . The SO(IO) model (in which the cancellation is compatible with quark/lepton mass 
hierarchy) provides the rough structure of proton decay operators given by YdYrf, and 
the size of and Cr operators are much larger than the current experimental bound. 
We also stress that in SO(IO) models, the Higgs triplets in 10 and 126 have opposite 
D-parity and therefore appear with opposite signs in the Cl,r expressions as in Eqs.Q 
and ( ITOi) . Therefore, the cancellation (between h and /) is unnatural. On top of this, 
such cancellation need to happen for each decay mode. 

2. Since the cancellation between h and / is not naturally realizable, it is preferable that 
h coupling structure is similar to the up-type quark mass hierarchy, and the down-type 
quark and charged lepton mass hierarchy is generated by / and h' matrices. Even in 
this case, the hijhki contribution is comparable to the current experimental bound even 
if hu ~ Uu and ~ Vc, where i/u and i/c are up and charm quark Yukawa couplings. 
Therefore to avoid unnatural cancellation which has to be implemented for each decay 
mode, the simplest choice appears to be to have hu ^ ?/„ and ^ Vc- This implies 
that the h coupling has no role in generating the 1st and 2nd generation masses and only 
3rd generation masses are generated by h; in other words, h coupling can be chosen to 
be a rank 1 matrix to leading order. 

3. If the h coupling matrix is rank 1, and 1st and 2nd generation masses are generated from 
/ and/or h' couplings, thus avoiding the need for cancellation between h and /. We 
however still need suppression of the contributions from / and h'. The simplest choice 
seems to be that the 1st generation masses and Cabibbo mixing are generated by h[2, 
and / coupling does not contribute to them. Since the h' matrix is anti-symmetric, it 
does not contribute to the kl part of the Cl operator Eq.(l9]). In fact, since the down 
quark Yukawa coupling is much too large to satisfy the current experimental bounds for 
the the nucleon decays, the question of how to generate the first generation masses and 
Cabibbo angle becomes an important one and this is the reason way we adopt a 120 
Higgs representation instead of an extra 10 Higgs field. 

4. The size of fu components are also important for a natural suppression of the proton decay 
amplitude. If the / coupling does not contribute to the down-quark and electron masses, 
and Cabibbo mixing angle, fu and /12 can be taken small. In fact, this choice of / and h' 



coupling matrices can be simply consistent to the empirical relations: sin^^; ~ ^rrid/ms, 
and mdnismb — rriem^mr- Under the assumption that the h coupling is dominantly large 
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and rank one and h' and / are the correction needed to generate the masses of the first 
and second generation quarks and leptons, respectively, the quark mixings Vcb and Vuh 
are generated to be small due to the left-right symmetry. 

5. The last piece in this discussion involves the magnitude of the 13 element of / coupling. 
In the case of type II seesaw (triplet part dominant), a sizable /13 is needed to obtain a 
proper solar mixing and neutrino mass ratio. In this case, one still needs a cancellation in 
the nucleon decay amplitudes. It is important to study whether the /13 coupling can be 
made to be small. The size of /13 is important not only for the nucleon decay suppression 
but also for the lepton flavor violation as we will see next section. 

4 Lepton Flavor Violation 

We now turn to the discussion of lepton flavor violating rare decays in this class of models. 
It is well known that the SUSY GUTs predict the lepton flavor violating rare decays due to 
the possibility of large flavor mixings in the superpartner sector. In order to avoid too much 
FCNCs induced by them, the flavor universality of the SUSY breaking mass parameters is often 
invoked. Even so, the loop correction via the heavy GUT particles can generate flavor violation. 
Though there can be new flavor changing effects in quark sector in SUSY GUTs, the recent 
LHC result imply no large effects [241 [25]. We, therefore, concentrate on the flavor violation in 
the lepton sector. 

The typical sources of LFV in SUSY SO(IO) follows: 

1. e'^u'^Hc coupling (which is lOm • lOm ■ H5 coupling in SU(5)). 

2. Dirac neutrino coupling (Y^iu'^Hu)- 

3. Majorana neutrino coupling (/^££Al). 

In SO(IO) model, the Dirac neutrino Yukawa coupling Y,y is the effective combination of 10, 126 
and 120 doublet coupling with appropriate mixing ratios (as in Eq. 6) whereas the Majorana 
neutrino coupling is unified to the 126 coupling (up to a Clebsch-Gordon coefficient). 

One can express the RGE-induced off-diagonal elements of SUSY breaking masses in terms 
of the Yu and Yukawa couplings as follows: 

A/Vf2 3mg + Al ^ 9mg + 3Ag ^ f,.)]r,MR (u) 

^g^^ ^(y,,,r.,,,)ln^^^ ^g^^ ^(/,,,/.,,,)ln^^, (11) 

The initial flavor universality is assumed and Mfj = mgl, and Aq stands for an universal 
trilinear scalar coupling. 
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In general, the off-diagonal elements by these couplings can be parametrized as 



k2 I U\ (12 ) 

V 1 



where the unitary matrix U is just for the parameterization in general. The mixing angles in 
U depends on the source of flavor violation. Once the source of the flavor violation is specified 
(with some assumptions), numerical quantities in the parameterization can be related to the ob- 
served values. For example, if the source only comes from the Dirac neutrino Yukawa coupling, 
= UlYI^^^^U^i (in the basis where charged-lepton mass matrix and the right-handed Majo- 
rana mass matrix are diagonal), we obtain U = Ul, n = {3 + ml/Al)/{8n'^){Y^"^^)l^ In Mu/Mr^. 
The quantities fcj are specified by the hierarchy of the eigenvalues of the Yukawa couplings (with 
RGE corrections). If Ur = 1 is supposed, the unitary matrix corresponds to the neutrino 
mixing matrix (up to renormalization group evolution). In general, the unitary matrix Ul is 
different from the neutrino mixing matrix. However, it is often assumed that the large neutrino 
mixing angles originates from Ul- In the SO(IO) model, the angles in U in the parameterization 
is different from the neutrino mixing angles, and those are computed from the expression of Y,^. 
In the triplet part dominant type II seesaw (i.e. type I seesaw part is negligible), the unitary 
matrix U corresponds to the neutrino mixing matrix if the fi, coupling is the dominant source 
of the flavor violation. 

The most stringent constraint is provided by the experimental upper limit on the branching 
ratio of — )• 67 [26] : 

Br(/i ^ 67) < 2.4 X IQ-^^ (13) 

To suppress the /i — >■ 67 amplitude in SUSY models, the 12 and 13 elements of the slepton 
mass matrices, given by the equation below, have to be small: 

= -/tm^ (^A;2 sin 2^12 cos 023 + e*^ sin 013 sin 023)e^^'^""\ (14) 

= -Kmg(^fc2 sin 2012 sin 023 -e^^ sin 013 cos 023)e'^, (15) 

where 6ij denotes mixing angles in the matrix U and a, /5, S are phases in U. 

The left-handed fi ^ e'-f amplitude can be cancelled by a choice of ^2 and mixing angles. 
However, such cancellation cannot happen for both yU — )■ 67 and fi-e conversion simultaneously. 
As a consequence, the suppression of /i — )■ 67 and /i-6 conversion requires (at least) one of the 
foUowings: 

1. small K (which means that the relevant Yukawa couplings are small). 
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2. both 12 and 13 mixings in U are small, 



3. ki 1 (which means that the eigenvalues of Yukawa matrix is hierarchical) and 13 mixing 
is small. 

If K is small, r — )■ fi{e)'j is also small. In the other cases ( (2) and (3) above), r — )• /X7 
can be sizable, staying just below the current experimental bound while satisfying the bounds 
from /i — 67 and fi-e conversion if 23 mixing angle in U of the source is large. It is therefore 
important to sort out the models to see if r — )■ /i(e)7 can be observed just below the current 
bounds, with a potential for /i — )■ 67 to be observed near future. We will study the patterns 
of LFV in our SO(IO) model using type I and type II seesaw models once we fit the neutrino 
masses and suppress proton decay. 



5 Predictions for flavor structures in type II seesaw 

We discuss the suitable choice of flavor structure in lepton sector, and how to accommodate 
two large neutrino mixing angles(for atmospheric and solar neutrino oscillations) and one (rel- 
atively) small 13 neutrino mixing angle naturally in type II seesaw in SO(IO) GUT |23j. 

Let us first describe the flavor structure in a general case. We start from a basis where 
neutrino mass matrix is diagonal. In this basis, the neutrino mixing matrix is equal to the 
diagonalization matrix of charged lepton mass matrix Mg. 

f/^MNsMeMlt/^MNS = {Mf^')'. (16) 

Since the muon mass is much smaller than the tau lepton mass, one can decompose as 

Me = M° + 5Me, (17) 

where is a rank 1 matrix and generate tau mass, and elements of SMf, is smaller than M^. 
The rank 1 matrix is written in general as 




In general we can make the rank 1 matrix symmetric by rotating the right-handed lepton 
fields, and therefore, we write it in the symmetric form. In that basis, 6Me is not necessarily a 
symmetric matrix. The rank 1 matrix can be diagonalized by two angles (or one can say one of 
the three angles is unphysical if 6Mf. = 0). A unitary matrix to diagonalize the rank 1 matrix 
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can be written as 



Uo 



/ 1 \ / cos^^ -sin^, \ 

cos^a — sin^a sin^s cos^g 
^0 sin 9a cos 9a / \ 1 y 

/ COS 9 s — sin 9 s 

cos 9a sin 9s cos 9a cos 9s — sin 
\ sin 9a sin 6'^ sin 9a cos 6'^ cos 9a 



where 



tan^s = -, tan^a = 
b 



(19) 
(20) 

(21) 



We define a basis rotated by Uq. (We attach "hat" to distinguish from the original one.) 

M,^M^ + 5Me, (22) 

where = diag(0,0,m3) {m^ = + 6^ + c^), and SMe = UQ{SMe)UQ . Define the diagonal- 
ization matrix of Me as Vg. Then the mixing unitary matrix is written as 



The 13 element of C/pmns is 



C^PMNS = VeUo- 



Ue3 = -(V;)i2sin6'„+ (V;)i3Cos6'a, 



(23) 



(24) 



Prom our assumption TO3 » {SMf.)ij, we expect (K)i3 and (14)23 to be small. The angle 9a 
is almost same as the 23 mixing (for atmospheric neutrino oscillations) if (14)23 is tiny. The 
solar angle is modified from 9s by (T4)i2- Naively, if (T4)i3 can be negligible, the 13 mixing is 
(K)i2/ sin^atm- At this stage, this is just a parametrization of 13 mixing. But it is very useful 
to work using the unification picture. Interestingly, the experimental measurements of the 13 
mixing are consistent with (14)12 = Vus, which is the Cabibbo angle in the quark sector. 

We note that C/e3 is exactly equal to zero in the limit where only 33 element of SMg is 
non-zero (in this limit, the rank of Mg is 2, and electron is massless). In general, 9a and 9s 
are large. If {(^M,^)^^ dominantly generate muon mass, the 13 mixing angle is naturally small 
compared to the others, and the naive size of the 13 mixing is expected to be electron/muon 
mass ratio. 

Generically, there is no reason why (5Me)33 is dominant compared to the other elements 
and (K)i2 is small (in the "hat" basis). In the triplet-dominant type II seesaw, however, such 
situation is quite natural because the correction of the charged Icpton mass matrix and the 



neutrino mass matrix can be unified (up to factor) to the 126 Higgs coupling /. Besides, the 
126 coupling is good to generate the strange quark and muon masses for the Clebsh-Gordon 
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coefficient. We emphasize that the situation is consistent with the suppression of nucleon decay 
amphtudes by a flavor structure. As we have explained, we are studying the situation where the 
fermion masses are dominantly given by rank 1 matrix (h), which gives third generation charged- 
fermion masses, and the / and h' coupling matrices give the first and second generation masses, 
and fermion mixings. Under this assumption, qualitative structure of the fermion masses and 
mixings can be easily reproduced in the case of triplet-part dominant type II seesaw. In addition 
to the natural realization of two large and one small neutrino mixings, the quark mixings are 
small under this assumption due to left-right symmetry in SO(IO). 

Let us illustrate the feature of the flavor structure in triplet-dominant type II seesaw. We 
choose a basis where / matrix is diagonal, and parametrize 

/ be ac 

be b' ab \, / = /33 I A2 I . (25) 




h = 

\ ac ab a? 

In this basis, the neutrino mass matrix is (nearly) diagonal neglecting the type I seesaw term 
(which is assumed to be suppressed by heavy right-handed neutrinos). The neutrino mixing 
matrix is the diagonalization matrix of Y^. For simplicity, the h' contribution is neglected in 
this illustration. Then, in the limits of Ai, A2 — ?■ and /sa <^ /133, one obtains 

6^ + c 
tan2 6'atm = 5 — , tan 6*^01 = 7, sin 6*13 = 0. (26) 

The limit Ai,A2 — )■ corresponds to the massless electron limit and two zero eigenvalues of 
neutrinos. Therefore, it is naturally realized that two solar and atmospheric mixings are large, 
and small smOi^. The size of the 13 mixing is expected to arise from the electron/muon mass 
ratio, as well as the ratio of mass squared differences, Am^g/^'^is- order to fit the electron 
mass without fine-tuning, the component is useful. If the fine-tuning to fit the electron mass 
IS absent , the 13 neutrino mixing can be directly related to the electron/muon mass ratio, and 

The first possibility is that the 13 neutrino mixing only depends on electron/muon mass 



mass ratio. This can be constructed to make the 126 coupling matrix tri-bimaximal form 
by using a discrete fiavor symmetry [23J. In this case, the Am^2/^'^23 dependence can be 
dropped, and the predicted 13 mixing is 

iTfl 1 

sin 6113 ~ ./— sin6'atm ~ -KsSin6'atm ^ 0.05. (27) 
y 3 



^ If one allows a cancellation to fit the electron mass, the prediction is lost, and the 13 mixing can have 0(1) 
factor ambiguity. However, if one fit all the fermion masses and mixings without 120 Higgs coupling (only 10 
and 126 couple to fermions jlOl [27]). there is predictivity of 13 mixing due to the reduction of parameters [11) . 
In this case, the electron mass is obtained by a cancellation. The prediction of the 13 mixing is related to the 
Cabibbo mixing angle (roughly Vus/^f^), and the prediction can agree with the recent measured value. 
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This does not match the measured 13 mixing angle. We note that the tri-bimaximal form of 
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126 is not preferable for the nucleon decay suppression, and we discard this possibility even 
without using the ^13 mixing. 

The second possibihty is that the 13 mixing angle is related to both the electron/muon mass 
ratio and l^m\2/ ^m\^. The Ami2/Am|3 dependence of the 13 mixing angle is related to the 
/ coupling structure, which is related to the proton decay suppression, which we will see later. 
Indeed, in the rank 1 structure, the size of /n and /13 are important for both proton decay and 
LFV. The small /n component is favorable to suppress proton decay amplitude, and it allows 
prediction of the 13 neutrino mixing angle. 

Let us study how the 13 neutrino mixing angle is numerically restricted in the rank 1 h 
Yukawa structure: 



(28) 



Note that we can parametrize /12 = /21 = without loss of generality, by diagonalizing 1st- 
2nd block of the matrix /. As we have noted, — >■ is preferable to suppress proton decay 
amplitude and to obtain the empirical relation T4s — \ffndjrris. Roughly speaking, the 3rd 
generation masses are given by hz, and the 2nd generation masses are given by /22 with r2 
and Clebsch-Gordon coefficient. The elements gives the first generation masses as well as 
the Cabibbo mixing angle. The absence of /n provides natural understanding of the so-called 
Georgi-Jarskog relation (without a cancellation): mem^mT- ~ m^msmb by a choice of |ce| — 1. 
Other parameters can be chosen to be consistent with the other quark mixings and atmospheric 
and solar mixing angles. 

Prom /12 = 0, we obtain a simple relation among the eigenvalues of / and mixing angles: 

h 

k 
h 



13 



u sec^ 



13 



13 



we^*^) sec tan Q^^, t^-ii ^12 tan ^23, 
ue^^^) sec tan cot tan 6*23; 



(29) 
(30) 



where 
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cos ~ sin 
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(31) 
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(32) 



In triplet-part dominant type II seesaw (and in normal hierarchy), we have 
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Figure 1: Magnitude of sin ^13 depending on (5/7r (and also ^23, each choice of which gives one 
hne), which is restricted in a narrow band. 
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Am^g 



\h?-\f2? 
In the hmit m — )• 0, we obtain 

4 sin^ Qlj, CSC 2e\^ tan 6*^3 (sin Q^.^^ cos b + cot 2615^2 tan 61^3) 



(33) 



Am|2 
Amia 



1 - sin^ ^1^3 (2 + cot^ tan^ Q^j) - 2 sin^ ^^3 cos 5 cot ^5^2 tan ^^3 



(34) 



As a result, ^13 is bounded from below by Am^2/^^23- Even if m 7^ 0, ^13 has to be bounded 
because we obtain /1//3 — > u, /2//3 — > u in the limit of ^^13 — t- 0. In figure 1, we plot the 13 
mixing angle as a function of the CP phase 5. Since the 13 mixing angle is sensitive to the 
deviation from the 45 degree of 23 mixing, we vary the 23 mixing angles. It is not sensitive to 
the 12 mixing angle. It is interesting to note that 5 ~ vr is preferable if the 13 mixing is larger 
value (depending on 6*23). 

The 13 neutrino mixing angle can be modified by the diagonalization matrix of charged 
lepton masses, but the size correction is less than 0.05 radianj^l- In conclusion, tiny ^^i3(< 0.05) 
is not allowed for the rank-1 structure with type II seesaw, irrespective of the detail of fitting of 
fermion masses and mixings. This is in contrast with the (nearly) tri-bimaximal model, where 
^i3~l/3KsSin^23~0.05. 

The suppression is 11,12 elements of / coupling reduces the // contribution of the nucleon 
decay amplitudes drastically. Because the / matrix is generated the large neutrino mixings 
and the neutrino masses directly in the triplet dominant type II seesaw, the 13 element of / 



^ The correction depends on relative phase freedom. The maximal correction (^ 0.05 radians ) occurs when 
the relative phase is or tt. 
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cannot be zero. Therefore, it can still contribute to the nucleon decay amplitudes and we need 
a cancellation which can be achieved by choosing the colored Higgs mixings to suppress them. 

As we have noted, the LFV can be also generated from the / coupling if the doublet Higgs 
mixing angle is small and the original / coupling in the superpotential is large since (at least) 
S\J{2)l triplet is hghter than the GUT scale. In that case, however, — > 67 is generated and 
the doublet Higgs mixing at GUT scale is bounded. If the doublet Higgs mixing is 0(1), and 
/ ~ /, the size of Br(/x — > 67) is the same order of the current experimental bound. Since the 
size of 13 element is predicted in the type II seesaw model, the Br(T — > /ij) can be predicted, 
and it is about O(10~^°), which is below the sensitivity at the LHCb. In the type II seesaw 
neutrino, it is not possible to enhance Br(T — > /I'j) to be of the order of the current experimental 
bound while satisfying the bound of // — > 67 and /i-e conversion. 



6 Type I seesaw 

In the case of triplet-part dominant type II seesaw, the simple formula of the mixing angles are 
obtained, and the two large and one small neutrino mixing angles can be easily realized. On 
the other hand, for the type I seesaw case, the situation is more complicated since the seesaw 
formula includes an inverse of the / matrix. In this section, we investigate how the fermion 
masses and mixings are reproduced in type I seesaw with suppressed nucleon decays. 
Because the Dirac neutrino Yukawa coupling is given as 

Y^^h- 3r2/ + c,h', (35) 

the type I seesaw neutrino mass matrix is proportional to 

YJ-'Yj ^{h + c,h')f-\h - c,h') - 6r2h + 9r^/. (36) 

Because the up-type quarks are more hierarchical than down-type quarks, r2 is small in the 

current scheme, and the last two terms can be negligible. We, therefore, concentrate on the 

first term N = {h + c^h')f-^{h - c^h'). 

We denote h and h' matrices as 

(hi \ / ci -C2 \ 

h=\ h2 , h' =\ -ci C3 . (37) 

V h J V C2 -C3 / 

We study whether a hierarchical / matrix required by charged fermion fits and also proton 
decay suppression is consistent with the profile of observed neutrino mixings. For this purpose, 
we first obtain the / coupling by solving N = {h + Ci,h')f ^^{h — Cyh'). Denoting 

N = C/diag.(ni, 712, n^)U^, (38) 
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we obtain 



Defining Xj as 



we obtain 



/ = (/i - c,h')U*diag. ( — , — , —) U'^ih + c,h'). 



U^'ih + c.h'), 



\ rp \ rp \ rp 

f = — Xi xi H X2 X2 H X3 X3. 

ni n2 ris 




The neutrino mass matrix is expressed as 



(39) 



(40) 



(41) 



(42) 



where Vy, is a VEV of up-type Higgs field, and vr is a VEV of 126 Higgs field which breaks 
SO(IO) down to SU(5). (More precisely, this / coupling is original / coupling as a convention). 
As we have noted, we are neglecting terms from 9r|/ — 6r2h, because each component of N 
is 0(100) as we will see. The unitary matrix U is the neutrino mixing matrix (up to the 
diagonalization matrix of charged lepton mass matrix) . 

One can express / matrix using a general form oiU. To capture the essence of the discussion, 
we first use a tri-bimaximal form for U, and later on add the corrections due to non-zero ^13. 
Using the tri-bimaximal form 



U 



/ ^ J- \ 
]_ j_ ]_ 

]_ j_ j_ 



(43) 



we obtain 



xi = -j={c^{ci - C2) + 2hi, c^(2ci + C3) - /i2, -c,.(2c2 + C3) - h^), 



X2 = 



X3 



1 

73 
1 

72 



'Cu{-Cl + C2) + hi, Cj.(ci - C3) + /l2, Cu{-C2 + C3) + hs), 

;c^(ci + C2), -c,,C3 + h2, -c^C3 + ha). 



(44) 
(45) 
(46) 



Our purpose is to find a solution where hi,h2,fu are suppressed. The solution can be found 
by assuming 



Ci — C2 = {rii, 712 '^'^3, normal hierarchy), 
ci + C2 — (ni, 77-2 ^ 7i3, inverted hierarchy). 



(47) 

(48) 
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In the case of normal hierarchy, we need vn/n^ ~ lO^'^ GeV, and therefore, we expect 
77-3 ^ 100 — 1000. Assuming rii <^ n2,n3, for the case of normal hierarchy, we have n2/n^ — 



^"^soi/^"^atm ~ 0.15. By a choice of Ci = C2, the / matrix is obtained (assuming hi, /i2 ^ 
ci, C2 ^ C3 ^ /la to express in short) as 



6ni 



— I -2c^CiC3f^ clcl c,c,h |. (49) 



2c^ci/i3g c^c^hs hi 

We emphasize that the 11,12,13 elements of / are suppressed by ni/ris, which is precisely what 
is needed to suppress nucleon decay amplitudes. 

In the case of inverted hierarchy, we need |mi| ~ 1^21, and |mi| — |m2| = Am'^^^^j. / {2mi) . 



If m3 <^ mi, nil — y ^fnltm- So, rii ^ 100 — 1000 in this case. Assuming ^ ^1,^2, we 
can obtain that 11,12,13 elements of / are suppressed by n^/rti similar to the normal hierarchy 
case. For example, the 13 element of / for the choice of C2 = — ci and /ii = is 

2Cj,Ci(-/l3 - Ct,C3 + 2CJ.C2) -2CyCi{h2. + C^yCi + CyC^) 

/i3 7. \ 7. , loUj 

brii 6n2 

which is equal to —CuCih^/ni (for rii ~ 71,2), CyCih-^/'ini (for rii ~ — ^2), for example. 

In both cases, the fu elements can be suppressed naturally, and the suppression of fu is 
related to the hierarchy between the VEV of 126 and v^jvuy ~ 10^^ GeV. 

We note on the effect of modification from the tri-bimaximal case : (f/13 7^ 0). In the 
inverted hierarchy case, the condition is just changed to C2 cos 6^23 + Ci sin 6*23 = (as long 
as hi = 0). In the normal hierarchy case, we require (xi)i = to make fu — )■ 0. In that 
case, (x3)i ~ ci similarly to the tri-bimaximal case. In the tri-bimaximal case, (x2)i = is 
satisfied (for hi = 0). In the case of 613 7^ 0, (x2)i ~ ci sin 6^13 /(sin 6*12 sin 6^23). Therefore, 
the fu elements are suppressed by Uf^n2/n3 and rii/n^. Because of the approximate relation 



?T,2/?^3 — y Am'^^Y/ ^''T^atm ~ ^13) the size of the fu elements is not far different from the case of 
013 = 0. 

We solved the / matrix which can reproduce the mixing angles in U by Eq.( l39|) . As we have 
seen, for the solution with suppressed fu, Ci and C2 has to be related (|ci| = \c2\ for the tri- 
bimaximal mixing). If we put the mixing angles in U, the condition is given as tan 6*23 — |ci/c2|. 
Inversely speaking, if we start from a / matrix with suppressed fu to suppress nucleon decay 
naturally, the large atmospheric mixing angle implies |ci/c2| — 1. Interestingly, it naively 
implies a "post-diction" for quark mixing : V^b ~ Vusfns/mi,. 

We note on the fitting of the charged fermion mass and mixing. The size of c^cs has 
to be 0(0.1) as long as the muon mass comes from /22. To obtain the proper value of Vet, 
we need a cancellation between / and h'. It can also modify the atmospheric mixing from 
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the diagonalization of charged-lepton unless Cg ~ +1 (For example, at the SU(5)-like vacua, 



Apart from the detail fit, we have showed that the suppression of fu and hi^2 is possible to 
reproduce the neutrino mixings in type I seesaw. Although the detail of the charged fermion 
masses and mixings may depend on the threshold corrections (both GUT scale and weak scale) 
or any other possible higher order effects, the suppression of fu can generate an interesting 
feature of the type I seesaw. The fu elements is suppressed by rii/ris, as we have described. In 
this case, the / coupling can be a source of r — /i7 without enhancing yU — )■ 67 and r — > 67 if 
one of the SM decomposed representations in the 126 Higgs fields is lighter than the unification 
scale [18]. The feature of this solution is obtained because the fu components can be tiny to 
realize the neutrino masses, which is suitable to suppress proton decay amplitude. This is the 
main difference between type I and type II seesaw. In the case of triplet-part-dominant type 
II seesaw, it is impossible to enhance r — t- yU7 around its current experimental bound after 
satisfying the bounds of yU — )■ 67 and fi-e conversion because a size of /13 element is needed 
to generate the proper neutrino oscillation data. In both cases, /i — )■ 67 can be generated just 
below the current experimental bound via the Dirac neutrino Yukawa coupling, or left-handed 
Majorana neutrino couplings. If r — )■ /i7 is discovered soon, the structure in type I seesaw with 
suppressed nucleon decay is preferred. 

7 Predictions of /i — > 67 decay 

In this section, we show the predictions of the branching ratio of /i — )■ e7 decay in the SO(IO) 
model with proton decay suppression. 

As is expressed in Eq. (fT2l) . the size of the off-diagonal elements is specified by k, and the 
numerical quantity of k is specified by the size of the coupling matrices and the SUSY breaking 
parameters. For example, if the source is the Dirac neutrino Yukawa coupling, the numerical 
quantity is roughly estimated as 



where Y^^^^ is the 3rd eigenvalue of the Dirac Yukawa coupling matrix, M^, is a cutoff scale. 
Mm is the 3rd right-hand Majorana mass, mo is a SUSY breaking universal scalar mass, and 
Aq is a universal scalar trilinear coupling. 

In Figl2], we plot the branching ratio of yU — ?■ 67 decay for various scenarios. We choose the 
unified gaugino mass to be M1/2 = 700 GeV (corresponding to the experimental lower limit of 
rrig ~ 1.7 TeV for gluino mass) to satisfy the recent LHC experimental bounds, and tan/3 = 50. 





(51) 
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10-11 




mo(GeV) 

Figure 2: Branching ratio of /i — >■ 67 decay vs. the universal scalar mass ttiq. We choose 
Ml/2 = 700 GeV, and tan/3 = 50. The universal trilinear scalar coupling is chosen to satisfy 
rrih = 125 GeV. The current experimental bound is Br(/i — )■ 67) < 2.4 x 10~^^ [26] . It is 
expected that the decay can be observed if the branching ratio is larger than 10~^^ at the MEG 
experiment. The detail of the plots are given in the text. 

In the plot, the universal trilinear coupling Aq is chosen to be a value to make the lightest Higgs 
mass to be rrih = 125 GeV [28j. For example, for the universal scalar mass mo = 1000 GeV, we 
obtain Aq = —2000 GeV (In our sign convention, the positive value of M1/2 takes the trilinear 
scalar couplings to the negative direction by RGEs. Namely, for a boundary condition Aq = 0, 
the scalar trilinear couplings are negative at the weak scale in the convention). We note that 
the flavor violation is enhanced because a large magnitude of Aq is needed to obtain the Higgs 
mass. For larger mo, the left-right stop mixing becomes smaller (for fixed Aq). Therefore, in 
order to obtain rrih = 125 GeV for a larger value mo, a larger magnitude of Aq is needed, and the 
flavor violation is more enhanced. As a result, the behavior of the plots are different from the 
simple dependence of the slepton masses. Due to a large magnitude of Aq and a large left-right 
mixing, one of the stau masses becomes tachyonic for a small value of mo. As a result, mo has 
to be larger than 650 GeV for the plot. For the given boundary conditions, squark masses are 
heavier than 1.5 TeV. 

We choose the following four scenarios to plot: 

1. (Red lines) 

The flavor violating source is the neutrino Dirac Yukawa coupling, and the unitary matrix 
U in Eq. (fT2|) is the same as the neutrino mixing matrix. Namely, the Dirac neutrino mass 
matrix can be written as 

n = UYf'^^, (52) 
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in the basis where the right-handed neutrino mass matrix is diagonaL We choose M^3 = 
10^3 GeV (red sohd hne) and 3 x lO^^ Q^y (^^^^ dashed hne). 

2. (Blue sohd hne) 

The flavor violating source is the Majorana coupling of the left-handed lepton doublets 
and SI]{2)l triplets for the type II seesaw, and the triplet part provides the dominant 
contribution of the light neutrino masses (Then, the unitary matrix U is same as the 
neutrino mixing matrix). As we have explained, the coupling is the 126 Higgs coupling 
in the SO(IO) model, and the second generation fermion masses are generated by this 
coupling in the current setup for the proton decay suppression. There is a free parameter 
for the 10 and 126 Higgs mixing. We choose f^^"^^ = 0.1 for the original 126 Higgs 
coupling without multiplying the Higgs mixing. The numerical quantity k is propotional 
to /s'^^, and the branching ratio is roughly proportional to k^. Therefore, for f^^^^ < 0.05, 
the branching ratio is smaller than 10~^^. Because /ss provides the second generation 
fermions mass, and /23 generates the quark mixing Vcb, one can expect that the branching 
ratio can be more than 10~^^ in this scenario. 

3. (Blue dashed line) 

The flavor violating source is the Dirac neutrino Yukawa coupling, Eq.(l6]), in the SO(IO) 
models. Contrary to the SU(5) models, the Dirac neutrino Yukawa coupling is restricted 
in the current scenarios of the SO(IO) models. There is one free parameter c^, in the 
Dirac neutrino Yukawa coupling. We choose = 2 (which is the value of the SU(5)-like 
vacuum, which is preferable since one obtains Ce = — 1 and the Georgi-Jarskog relation can 
be naturally obtained). In order to observe the /i — )■ 67 decay at the MEG experiment [26], 
one needs a larger value of c^. The branching ratio is roughly proportional to c^. 

In addition to the above sources, there can be a flavor violating source in the e'^u'^Hc cou- 
pling, and the colored Higgs loop can induce the off-diagonal elements in the right-handed 
charged slepton mass matrix, and can generate LFV via neutralino loop diagram. The contri- 
bution is calculated to be less than 10~^^ for the branching ratio using the boundary condition 
of the SUSY particle spectrum. 

In the above calculations, we assume the gaugino mass unification and the universality 
of the scalar masses. However, if there is a SUSY breaking contribution from the anomaly 
mediation, the gaugino mass unification can be relaxed even in the GUT modes. In that case, 
the gluino mass bounds from the LHC experiments does not necessarily restrict the Wino and 
Bino masses, and the Br(/i — > 67) can be enhanced. 
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In summary, in the SU(5) GUTs with type I seesaw, the branching ratio of the // — > 67 
decay can be as large as the current experimental bounds, and its magnitude depends on the 
right-handed neutrino Majorana mass. In the type II seesaw SO(IO) model with proton decay 
suppression, it is expected that the — > 67 is observed at the MEG experiments. In the type I 
seesaw SO(IO), on the other hand, the branching ratio is smaller than 10~^^ (for a natural size 
of Ci,) and // — > 67 decay may be difficult to be observed. 

As is explained, the t ^ /I'j decay width cannot be large enough to be observed near future 
in type II seesaw SO(IO) (and type I seesaw SU(5)) models, satisfying the experimental bounds 
of // — > 67 and /i-e conversion. In type I seesaw SO(IO) with proton decay suppression, the 
T ^ /I'j decay width can become large by the 126 Higgs coupling without enhancing /i — > 67. 
The T — > 67 decay is not enhanced in both type I and II seesaw scenarios. These features are 
important to distinguish the GUT models and the vacua of the GUT symmetry breaking. 



8 Partial lifetime of nucleon in type I and type II 



The hierarchical structure of the Yukawa coupling matrix is needed to suppress proton decay. 
As we have studied, natural suppression of the proton decay amplitude is possible in the type 
I seesaw relating to the hierarchy between the GUT scale and the typical seesaw scale ~ 10^^ 
GeV. In the type I seesaw, the structure is really simple and it allows predictions for the partial 
decay amplitude, which we investigate in this section. 

We denote the coupling matrices in the basis where h is diagonal: 
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I hi 
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(53) 



As we have explained, we assume that the h coupling is rank 1 (namely, hi and /12 are irrelevant 
to fit fermion masses). Then, we can choose ^; = without loss of generality. In that basis, 
— > is important to suppress proton decay amplitude. In type II seesaw, x cannot be small 
to fit the large solar neutrino mixing. In type I seesaw, on the other hand, x is small and 
the smallness of x gives a predictivity to the decay amplitudes. In the following, we take 
hi,h2,u,v — )■ 0, but we keep x to describe the difference between type I and II seesaw. 

The left-handed proton decay amplitudes (from the LLLL dimension- five operator Cl) from 
the chargino dressing diagram can be written as 



Al ip Kvr) ~ -^pqI {xpi - C2P^)y cos Oc sin OcXp^Kp, 
Al{p KVf,) ~ -l^pglcip^y cos Oc sin OcXp^^Ku, 



(54) 
(55) 
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where /3p is a hadron matrix element of proton, 6c is the Cabibbo angle, pi and ps are the 
coefficients from the colored Higgs mixing for // and h'f contributions, respectively: 

XLl = XaA^Ya^ = {MT%i, (56) 

= V2Xaz^Ya^ = V2{M^%^. (57) 

Here, X and Y are the diagonalization unitary matrix of the colored Higgs mass matrix Mt 
and we used the same notation in Ref.|T7]. The factor Xp^^u includes a loop function and 
chargino mixing angles 6'„, 9^: 

Xp^Ku = {Ai - A2)(cos2 e^Hl^ + sin^ O^Hl^) + [Ai + A2)(cose„ cosO.Hl^ + sin0, smO^Hl^), 

(58) 

A, = 1 + ^{F + \d), A, = ^\d, (59) 

where F ~ 0.48 and D ~ 0.76 are coupling constants for interaction between the baryons and 
mesons, mp is the proton mass, and ms' is an averaged baryon mass ms' ~ ms ~ m\. The 
loop function is defined as 

where 

„i:.,y)^^(-J^-«J^], (61) 
X — y \x — 1 y — I J 

Hud is defined similarly by replacing the slepton mass into down-type squark mass, and m^^ is an 
eigenvalues of chagino masses. We assume that SUSY breaking squark, slepton mass matrices 
are proportional to identity matrix for simplicity. We neglect the subleading contribution from 
Vcb and Vub- The p — )■ Kug is suppressed by a factor ~ 9c, but can be generated due to the 
mixing between electron and muon in the basis where h is diagonal. 

We find that A^^p — )■ Kut-) cx xpi — C2P5 and Ai{p — > Ku^) oc Cip^. Therefore, if x — )■ (as 
in type I seesaw), the left-handed amphtude of p — )■ Ku can be suppressed by choosing a small 
P5. The smallness of p^ is related to a vacuum selection of SO(IO) breaking vacua. In type II, 
on the other hand, x cannot be small, and one needs to choose small values both pi and p^ to 
suppress proton decay. If there are numbers of parameter in the colored Higgs mass matrices, 
the suppression is possible. However, it is not related to the vacuum selection. 

Since in the case of x — )■ in type I seesaw case the cancellation between the Yukawa 
couplings for different Higgs representations (// and h'f for example) is not required, the ratio 
of partial lifetime is predictable. In fact, the decay amplitudes of to anti-muon are obtained as 

^L(p->7rV^) ^ (3pglciP5y sin Oc sin OucAsXp^^+x, (62) 
AL{p-^K^fi^) ~ -(3pglciP5y cos 9c sin eucA4Xp_,^+x, (63) 
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where 

= cos ^„ cos e,{Hl + //^J + sin 9^ sin ^.(i/i + i/^J , (64) 

and 

1 Tfl 

^3= + F + A4 = 1 + ^(F-L>). (65) 

V 2 TnB' 

The mixing angle is an angle between -u, and c quarks in the basis and therefore, 9uc ~ 
^Jmu/rric- Therefore, we obtain the ratio of partial decay width in type I seesaw as 

(66) 

The maximal atmospheric mixing in the type I seesaw requires |ci| ^ |c2|, and the ratio is 
predictive (up to the sfcrmion mass spectrum). In the case of type II seesaw, there is an 
additional parameter (a; — C2P5/P1), and the ratio cannot be predicted. We comment that 

r(p^xv+) " Al ^^^^ 

is obtained for both type I and II since only h'f contribution is dominated. For the decay to 
the third generation lepton (because tau lepton is heavier than proton, the decay to the Vr is 
the only case), // piece also contributes if /13 — x is not suppressed. Therefore, the ratio of 
partial decay width to neutrino and anti- lepton is not predictable in the type II, while it is 
predictable for type I seesaw. The same situation occurs for n — > ttP decay. 

In the above expressions, we neglect the contribution from the right-handed dimension- 
five proton decay operator Cr. They can contribute to the decays into P,-, and can alter the 
prediction of the ratio of partial decay widths. The Cr contribution (for hi,h2,u,v,x ^ 0) can 
be obtained as 

Ar{p Kur) = apycyr{AiCiC2XRs cos 6c + A2{ciC3Pr8 - C2ypR7 + C1ZPR5) sin 9c) 

X (i^^g cos 9u cos 9y + iJ^g sin 9u sin 9^) , (68) 

where ctp is a hadron matrix element, and 

PR5 = -V2Xa2^{Y,,-V2Yae)^-V2iM^'),2 + 2{M^\2, (69) 

PR7 = -\/2X„4T^(F„3-i;2) = -\/2((M^l)34-(M^l)24), (70) 
PR8 = 2X,2^{YaS-Ya2)^{MT%2-{M^%2. (71) 

For a large tan/? ~ 50, since q and y, z have to be large to fit the down-type quark Yukawa 
couplings, Cr contribution can make the decay width comparable to the current experimental 
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bounds, and the simple relation shown previously can be disturbed by the second term {A2 
contribution). If the Higgsino mass is much heavier than wino mass, the Cr contribution is 
suppressed rather than contribution. We note that if we do not assume hi, h2, u,v ^ 0, the 
magnitudes of Cr contribution can excess the experimental bounds even for tan ^ ~ 5 and we 
should not neglect the Cr contribution. 

9 Conclusion 

We study the Yukawa texture in a renormahzable SUSY SO (10) GUT model for neutrino 
masses that gives a "large" ^13 while at the time suppressing nucleon decay without invoking 
cancellation between coupling parameters. We consider cases with both type I and type II 
seesaw separately. In the type II seesaw scenario, we find that ^13 > 0.05 radians irrespective 
of the detail of fitting of fermion masses and mixings (i.e. independent of the number of 
parameters) and that the measured value of 13 mixing angle prefers 5 ~ tt for the CP violating 
phase in the neutrino oscillations. We then study the predictions of these scenarios for lepton 
fiavor violation. The branching ratio of — > 67 can be as large as the current experimental 
bound, and should be observed very soon. The Br(T — > /I'j) is predictable, and is about 
O(10~^°) in the type I seesaw case, if the bounds of // — > 67 and /x-e conversion are satisfied. 

In type I seesaw scenario, observed ^13 can be accommodated naturally. We find that 1-2 
and 1-3 fiavor violation from the FCNC source is suppressed compared to the type II scenario if 
the proton decay is suppressed naturally. This gives a larger value of Br(r — > ii'y) while — > 67 
and T — > 67 are suppressed. If the t ^ fj.'j decay is observed soon then definitely it will point 
towards a type I scenario in this kind of SO(IO) model and the prediction for Br(// — > 67) can 
be checked for confirmation. Detailed numerical study of the model for fermion mass fits as 
well as lepton fiavor violation is currently under way. The ratio of the partial decay width of 
proton to K/i and KP is predicted (up to uncertainties from sfermion masses) in this scenario. 

The next generation of the baryon number violating nucleon decay experiments at Hyper- 
Kamiokande along with the information about the low energy SUSY states from the LHC may 
be used to distinguish between the different scenarios. 
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